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In the study of QCD dynamics, C* boundary conditions are physically relevant in certain 
cases. In this paper we study the implementation of these boundary conditions in the 
lattice formulation of full QCD with staggered fermions. In particular, we show that the 
usual even-odd partition trick to avoid the redoubling of the fermion matrix is still valid 
in this case. We give an explicit implementation of these boundary conditions for the 
Hybrid Monte Carlo algorithm. 

1. Introduction 

Lattice QCD simulations are usually performed with periodic boundary conditions 
(BC). However, other type of BC may be important in certain cases. For example, 
a comparison between systems with different BC can be used to understand finite 
size-effects in lattice QCD. Some years ago 1 C-periodic BC were studied as an 
alternative to periodic conditions. Then they were considered with the general idea 
of studying the spontaneous symmetry breaking aspects of the QCD dynamics in a 
simple way. 2 In that work an analysis was done in the continuum, and it was shown 
that in pure gauge theory these BC break the Z(3) symmetry explicitely, which has 
important consequences for the high-temperature deconfinement phase transition. 
These conditions are also useful in numerical lattice simulations of this transition. 
When quarks are present, C-periodic BC break both chiral and flavour symmetries. 

These boundary conditions are also especially important when topological prop- 
erties are relevant in the system under consideration. This is the case of the lat- 
tice studies of confinement through monopole condensation. Recently the role of 
monopoles in connection with colour confinement has been evidentiated in SU(2) 
and SU(3) gluodynamics, 3 for which a disorder parameter based on the magnetic 
U(l) symmetry has been constructed and studied by Monte Carlo techniques. The 
disorder parameter is the vacuum expectation value (vev) of a disorder operator, 
which is an operator that creates a magnetic monopole in the gauge configuration. 
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The definition of this disorder operator requires C-periodic BC in the time direc- 
tion. For the pure gauge case, this means that the links at time t + N t , where N t 
is the temporal extension of the lattice, are the complex conjugate of the links at 
time t. The effect on the simulation algorithm is a simple redefinition of the staples 
containing links that pierce the temporal boundary. The natural extension of the 
procedure used for the pure gauge case to full QCD requires the implementation 
of C-periodic BC in the presence of fermions. In particular, we will be concerned 
with the case of staggered fermions. C-periodic BC modify the fermionic matrix, 
and many proofs of properties used for the setup of standard simulation algorithms 
no longer hold. 

C-periodic BC in the continuum are defined by the action of the charge conjuga- 
tion operator C on the fields. 2 However, lattice fermions are different from fermions 
in the continuum. In particular, C is not a symmetry of the lattice action with 
staggered fermions. It also breaks translation invariance for finite lattice spacing. 
However, there is a discrete symmetry of the staggered fermion action, 



(1.1) 



(here i indicates the lattice point, [7j i(U is the SU(3) matrix associated with the 
link leaving the i-th lattice point in the fx direction, tpi is the staggered fermion 
field at the point i, and r\i^ the usual staggered fermion phase), which corresponds 
to charge conjugation in the gluon sector but which in the continuum limit also 
contains a flavour transformation. 2 ' 4 We will call C* this symmetry of the lattice 
action: 

C *U^ = U*^ c *i/> i = e i $?, c *$i = -il>Je i , (1.2) 

where ej = (— \Y i+Vi+Zi+u , (xi,yi, Zi,U) being the lattice coordinates of point i, 
the "T" represents the traspose operation, and we will use ip* = V> T . Translation 
invariance implies that a BC must correspond to a symmetry of the action. C*-BC 
are defined as the boundary conditions corresponding to the symmetry (1.2): 

$ l+Ar - c *$ 4 , (1.3) 

where $ is a field, or tp, and N is the number of lattice points in the direction 
in which we use this boundary condition. 

In the chiral limit the lattice action has a J/(1)e <£> U(l)o chiral symmetry of 
independent rotations on (x + y + z + i)-even and -odd lattice points. C*-BC break 
cxplicitcly this symmetry to L/(1)e=o*- Baryon number J7(1)e=o is also broken 
cxplicitcly down to Z(2)e=o- 2 

In this paper we will have in mind the physical problem mentioned above: com- 
putation of the vev of a monopole creation operator in lattice QCD. This means 
that we will consider imposing C*-BC in the time direction, and periodic BC in 
the spatial directions. However this is just to fix the notation in what follows; the 
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C* conditions could in fact be assumed in any direction. The purpose of the paper 
is to show the theoretical framework to be used in a lattice simulation with these 
BC (section 2) and how usual algorithms 5 need to be modified (section 3). Our 
conclusions are summarized in section 4. 

2. Mathematical description 

Let us consider the partition function of lattice QCD with staggered fermions 

Z = J(VU)(V^V^)e- s ^ s f , (2.1) 

where S g (U) is the Wilson action for the pure gauge sector, and Sf is given by 
Eq. (1.1). The fcrmionic variables can be integrated out to give 

Z = j (VU)e- s ^ u) det M(U) , (2.2) 



where* 



M{U)^ = mSij +J2l (EMj-m - Ul^Ju+J (2.3) 



2 



is the fermionic matrix, and periodic BC are assumed. Using this matrix notation, 
we can write 

S f = 4>Mi}j = i (ipMip) + {^Miff , (2.4) 

since Sf is a number. We can use the new variable defined as the column vector 
formed by ip and ip* , so that the fermionic integral is 

J (Zty>)(£>V0 e- s f = J (Dtf) e -i* T ^* = Pf(A) , (2.5) 
* T A*^( ^ j, )a(^ ) , (2.6) 



> i ,v T i • (2.7) 



with 

A = ( 

x M 

and we have introduced the Pfafnan of the matrix A, Pf (A). It is well known that 6 

Pf 2 (A)=detA. (2.8) 

Now let us consider C*-BC in the time direction. Then, following Eqs. (1.2) 
and (1.3), ipjy t = eoV'o; an d V'jVt — — e o'</'o\ where we have written in the subscript 



*Wc have absorbed the staggered phases by a redefinition of the link matrices U. 
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the temporal coordinate and omitted the spatial coordinates. Now Eq. (2.4) gives, 
for the terms connecting the slices N t — 1 and N t , 



2 (V 5 



2 

\ (V^ t _ M V^ t _i + 4>N t -iUN t -i,t4o) eo • (2.9) 



In this way the matrix A of Eq. (2.7) is substituted by 

A =( B M -.%)• 

where B satisfies the properties: 

B^ = -B* B = -B T (2.11) 

and M is the fermionic matrix Eq. (2.3) apart from the terms connecting the slices 
N t — 1 and N t , which have gone to the matrices B and — B*. 

Eq. (2.5) is still valid and we are interested in calculating Pf(A). 

2.1. Pseudofermionic variables 

The usual approach 7 to the simulation of theories with dynamical fermions is to 
rewrite the determinant of the fermionic matrix in Eq. (2.2) using that 

det(M t M)cx J V<fJV(f> expI-^MtM)- 1 ^] , (2.12) 

where is a complex bosonic field with the same quantum numbers as the Grass- 
mann field. One introduces the matrix (M^M) -1 , instead of M^ 1 , so that the 
pseudofermionic fields can be generated using a simple heatbath method. Since 
detM is a real number, det(M^M) = (detM) 2 . So, actually, this corresponds to a 
double number of flavours with respect to the original theory. However, the matrix 
M'M has two important properties: it has no matrix elements connecting even 
and odd lattice sites, and the determinants of its submatrices on the even and odd 
sites are equal. 8 Therefore one can avoid the redoubling of flavours by defining the 
pseudofermionic field only on even lattice sites. 

A remarkable difference between the partition function with periodic BC and 
the partition function with C*-BC is that in the latter case the determinant of 
M, Eq. (2.2), has to be replaced by the Pfaffian of A, i.e. ±\/dct A. Because of 
the square root, the usual trick of introducing pseudofermionic fields and rewriting 
this factor as the integral over these fields can only be applied if the number of 
continuum fermion flavours is such that the square root cancels. Moreover, in order 
to have a positive-definite integration measure, we need that the sign in front of this 
factor be +. Both conditions are satisfied if the number of continuum flavours is a 
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multiple of eight. This generates a further unavoidable redoubling. Until Sect. 3.2, 
we will be concerned with the numerical simulation of dot A, which is equivalent 
to simulating a double number of fermion flavours in the continuum limit (eight 
instead of four). In Sect. 3.2 we will discuss how to deal with the usual case of four 
staggered fermion flavours. 

Since the case of a system with eight fermion flavours in the continuum limit 
and C*-BC and the case of a system with four fermion flavours in the continuum 
limit and periodic BC are similar, we would like to follow in the former case the 
standard procedure to obtain the determinant of the matrix (2.10). First, in the 
Appendix A it is shown that det A is a real number. So we can also in this case use 
the matrix A^A to introduce the pseudofermionic field, which will have now twice 
the number of components as in the usual case. Second, we will now see that the 
matrix A^A does not connect even and odd lattice sites. 

Using the form of the fermionic matrix (2.3) and the definitions of B and M 
given by Eqs. (2.9) and (2.10), we can split the blocks of the matrix A in even and 
odd lattice sites: 



A = 



( \B eo -m -\Dj e \ 

\B oc -\Dj Q -m 

m \D co -\B* eo 

V \D oc m -\B* oc J 



(2.13) 



From Eqs. (2.3) and (2.11) it is easy to see that 



Dl = -D co Dl = -D* co (2.14) 
Bl = -B* eo Bl - -B co . (2.15) 

In the expression (2.13), the new blocks divide the matrix in rows and columns 
identified by the pairs [0,e], [0,o], [l,e] and [l,o], where the number indicates the 
blocks defined in Eq. (2.10) and the letter the even/odd subblock. Reorganizing the 
rows and columns to [0,e], [l,e], [0,o] and [l,o] (the determinant does not change), 
we rewrite the matrix A in the following form: 



where 





±4 

2 co 


A c 


M 





— m \ 


m 


o ) 


B co 




D co 


~B* eo 



(2.16) 



(2.17) 



(2.18) 



and A oe has the same form as A eo with the exchange e <-> o. Eqs. (2.14) and (2.15) 
give 

Al - -A* co . (2.19) 
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From Eqs. (2.3) and (2.9) we obtain the explicit form of the matrices D co and 
B co : 

(D co )j. k = ^2{Uj,»S kJ+ll - Uj^^Skj-p) , (2.20) 

where j is an even site and the link between j and k does not connect the t = N t — 1 
and t = time slices (otherwise (f C o)j,fe = 0), and 

(Beokfc = C/ fe f , t , ? = 0,fc = 7V t -l ^ (2.21) 
(B co )j.k = in every other case J 

where j = Nt — 1 or j = means that the temporal coordinate of site j is Nt — 1 
or zero, respectively. When j is an odd site, (D oe )j : k has the same expression as 
Eq. (2.20), and 

(Soekfc = -Ult j = N t -l,k = ] 

(Boehk = -Ut tt j = 0,k = N t -l \ (2.22) 

(B oe )j_k = in every other case J 

These expressions completely determine every element of the matrices A eo and A oe . 
We now compute 

A i A = / -j? + \Al e A oe -I(^ eo -4 eM )\ 

But using the properties (2.14) and (2.15) it is direct to see that the blocks out of 
the diagonal in Eq. (2.23) are zero. Then, using Eq. (2.19), 

and, as a result, we see that the matrix A^A connects only lattice points of the 
same parity. 

In order to use the same trick to avoid the flavour doubling produced by the 
introduction of A^A, that is, to define the pseudofermionic field on even sites only, 
we need to show that, also in this case, the determinant of the even and odd parts 
in the matrix (2.24) are equal. This is done in the following subsection. 

2.2. Reducing the number of flavours: even-odd partitioning 

Let us write K = A^A. Then we have 

_ / m 2 / 2 - \A* co A oc \ / K c \ 

K -{ m 2 / 2 - \A* oe A eo ) - { K a J ' (2 - 25) 

where I2 is the 2x2 identity matrix. We need to show that detif e = detif . The 
strategy will be to show that both det K e and det K a are equal to det A. To this 
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aim, let us consider the matrix A written in the e-o partitioned form, Eq. (2.16), 

(2.26) 

where 

(2.27) 

Ko=\ n°° 7." )■ (2-28) 



±4 


2^00 
(J, 


( ° 


—m \ 


\ m 


o J 


( B co 
V D co 


D* co 


( B oc 


D* 

oe 



and 

a ~'{d2 -$.)■ <2 - 29) 

We will make use of a general property of the determinant of a square 2N x 2N 
block matrix, 

det ( W ^Z^j = dct W dct ( Z ~ wx ~ ly ) . (2-30) 

where X, Y, W 7 Z are N x N square matrices and X is invertible. Applying this 
equality to the computation of det A, we obtain 

det A = det(fi) dot (V - j4>eM _1 A e ^j = det (V - ^lA^^A^j . (2.31) 
Using the explicit form of \i and A oc it is easily shown that 

Moc = -A* eM • (2.32) 
Inserting this equality in Eq. (2.31), it follows that 

det A = det (^i 2 + ^A* oc A co ^j = dot K , (2.33) 

where the last equality is implied by the fact that N is an even number. 

We can compute again det A using the same procedure after exchanging even 
and odd variables in A. Using a property analogous to that in Eq. (2.32), namely 
nA eo = —A* eo n, it is then easily shown that 

det A = dot U 2 + -A* co A oc ) = det K c . (2.34) 



4 

This is a proof that det K c — det K . 

We have shown that det K e det K a = (det A) 2 . On the other hand it is also true 
that det K e det K = dct(A' < A) = det det A, therefore we have also obtained a 
proof that det A = det A^ alternative to that given in Appendix A. Moreover, we 
have that 

det A = det K c = dot ( m 2 I 2 + -A f oc A oc > 0, (2.35) 
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since both m 2 / 2 and Al e A oc are positive definite matrices. Therefore detA is a 
positive number. 

We notice that this method can be easily applied also to the standard case, thus 
providing a proof alternative to those presented in Ref. 8. This is shown in detail 
in Appendix B. 

3. Hybrid Monte Carlo Implementation 

We will show now how the standard Hybrid Monte Carlo (HMC) algorithm 5 needs 
to be modified to incorporate C*-BC. In this algorithm one introduces fictitious 
momenta, conjugate variables of the links, as dynamical variables, and makes fields 
evolve with a mixed dynamics, in which deterministic and stochastic steps are alter- 
nated in a prescribed way. In the deterministic part of the algorithm, the system fol- 
lows the equations of motion derived from the Hamiltonian of the (4+l)-dimensional 
system. These equations give the evolution of the fields in the fictitious time r. The 
equation of motion for the matrix Uj tfl is 

where the conjugate momentum Hj^ is a traceless Hermitian matrix, and U is the 
derivative of U with respect to r. The equations of motion for the momenta H are 
obtained by imposing that the Hamiltonian be constant. The integration of these 
equations of motion is carried out numerically after discretization of r. Usually the 
temporal step is of order 10 -2 -10 -3 and the configuration space is sampled with a 
total length of the trajectory of order 1. The stochastic part of the algorithm consists 
in the generation of new momenta and new pseudofermionic variables according to 
their probability distributions at the beginning of each trajectory. Moreover, at the 
end of each trajectory a Metropolis accept-reject step is performed, which makes 
the algorithm exact. 



3.1. HMC algorithm with C* boundary conditions 

Once introduced the pseudofermionic fields, defined only on even sites, and the 
auxiliary momenta fields, the partition function of the system is 

Z = j (VUV^V^VH) e~ n , (3.2) 



with 



H = \ E trH L + S 9 + ^(^A)-^ . (3.3) 



2 

3,1* 



To obtain an equation of motion for H we require that TL be a constant of motion, 
that is, H — 0. The tricky part in this differentiation is in the fermionic term. The 
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derivative of Hf = ft (A* A)' 1 (j> is 



H f = -^(^(A^A)- 1 : 



dA ^ H JJ A A- dA * f/ T /T T A + 



dU. 



4t 9A u JJ. I /it dA TjT rjT _ rrT rj* dA a 



3,V 3,V fjTJ* 
3,H 



x(A f A)-\ 



(3.4) 



where we have made use of Eq. (3.1) and the fact that H is Hermitian. We have 
also used that A is linear in U, U T , U* and U\ and then the partial derivatives in 
the previous expression commute with H and U. 
It is convenient to introduce the operator 



P-ij — XiXj 



where 



X = . 
It is easily seen that P is Hermitian: 

(P% = (Pjrf = XiX] = Pij , 



(3.5) 
(3.6) 

(3.7) 



and, being <j> defined only on even sites, is taken to be zero unless i and j are 
both even sites. On the other hand, we have that 

/c^ty _ dA ( dA V _ 9A1 
( dAW _ dA ( dA \ f _ dA* 



(3.8) 



and then we can write 



3,H 



iHi 



jj. AP^ + U PA^ 



+ 



^IJap^u^ + pa^u^ 

\ 3,f 3,1* I 



H.c. 



(3.9) 



where H.c. means the Hermitian conjugate and we have used the cyclic property of 
the trace operation. 
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The trace in Eq. (3.9) is taken over both color and site indices. Taking this into 
account and using the cyclic property again, we have 



Hf = -£tr (H. 



3,V 



AP- 



0A1 



AP 



dUi 



PA f 



Ui 



PA"* 



dA 



dUi 



+ 



dA 



H.c. 



(3.10) 



The following step is to calculate the derivatives appearing in Eq. (3.10). From 
Eqs. (2.16) and (2.19), we have that 

At = -A* . (3.11) 

We write again A and A' in terms of subblocks, identifying the rows and columns 



A = 



, [0, o] and [l,o]: 






/ —TO 


1 R 

2 -Deo 




TO 


2 -^co 


1 D* 

2 co 


Id lm 

2 -Doc 2 oe 





— TO 


\ 2^°° — 2^oc 


m 





At 



V 





-rr. 
-D* 

2 oc 



— TO 
1 



m 


Id 

2 ^oc 



1 r>* 

— -D* 

2 eo 


— TO 



2 -°co 
TO 

o / 



(3.12) 



(3.13) 



Once written everything in terms of B and D we can easily see where U and U T 
appear, from Eqs. (2.20)-(2.22), and compute the derivatives. The result is given 
in Table 1. 

Now it is easy to calculate the different terms in Eq. (3.10). Let us do, as an 
example, the first one. The trace over lattice site indices affects only the expressions 
between parentheses. From Table 1, we see that, for j even, and j ^ N t — 1 or /i ^ t, 
the first one gives 



tr AP 



gAt 
9UZ 



= A 



k^kXkX^ (afj—) - ~9 Al,ob-+M,fc X fc X [0,c]i - 

\ OU 3^/k u k 2 Z 



-no _ 



--^(D oe )j + ^k X k + - (B* e ) j+ ^ k Xlj (Xt) 

-\{{A oe)j+ ^ k X k )\X^ = -\{A oe X)) +lx {X^ , (3.14) 

where we have used Eqs. (3.12) and (2.18), and called X® the first three components 
of the vector X k , and X^ the second three components. Notice that there is not 
a mass term because it is diagonal, which would imply that k = j + /i, that is, k 
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Table 1. Value of the derivatives appearing in Eq. (3.10). fi = t indicates the time direction, in 
which we are imposing the C**-BC. 



j even 



j odd 



j¥=N t -l 

or fl t 



j = N t -l 
and fi = t 



■\ ^a,[l,e]j 5/3,[l,o]j+ M 



I <W,o]j <*/3,[l,e]j+,u 



(ML) 



or fl t 



~\ Sa,[0,e]j S/3,ll,o]j+n 



-\S a ,[0,o}j$(3,{l,e]j+H 



j = N t 
and fi - 



-\ 8a,{0,e]j 8/3,{0,o}j+fi 



dA 
dUT 



i<P / a,f3 



or fi ^ t 



\ 5a,[0,o]j+ M ^/3,[l,e]j 



-| 4*,[0,e]j+/x <W,o]j 



and /i : 



\ 5q,[1,o]j+/x <W,e]j 



| 5 a ,[l, e ]j+/i ^,[l,o]j 



or fi t 



\ S a ,[l,o]j+fi 8/3,[0,e]j 



j = N t 
and = 



| <^q,[0,o]j+/x <^/3,[0,e]j 



-| ^a,[0,c]i+ M ^/3,[0,o]j 



would be odd, and then Xk = 0. For j even and j = N t — 1, /i = i, an analogous 
calculation gives 

This term does not contribute when j is odd, because in this case {X^)® is zero (X 
is defined on even sites only). 

Computing every other contribution in Eq. (3.10), the final result can be written 

= - E tr { H i# ( iF ^ iF L) } > ( 3 - 16 ) 

i,n 

where Fj^ takes a different form depending on j and fi. 
For j even and (j, fi) ^ (JV t - 1, t), 

Fi* = Uj,, (-\{A oc X)) +ll {X^ (±X]({A oc X)^ + ^ ; (3.17) 
for j even and (j, fi) — (N t — 

Fi* = ^> (-J(^oeX)^(Xt)0 + Qxj((A oe X)t)} + ^^ ; (3.18) 
for j odd and (j, n)^(N t -l,t), 

F^ = U jt , Ux^((A oc X)^+ (^(A oc X)° 3 (X^ +f }j j ; (3.19) 
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for j odd and (j, u) = (N t - 1, t) 



-X] + ^{A ocX y)) ( -{A oc X)%X^) 



(3.20) 



Since we determined completely the elements of the matrices A co and A oc in sec- 
tion 2.1, these expressions allow us to obtain as a function of the links U for 
every j and fi. 

Coming back to Eq. (3.3), the condition H = gives 

w = = E tr ( / ///- // /./- + § (m^u^v^ + h.c.) - <///,,/•;,, + h.c.; 

(3.21) 

where V} ljU is the sum of staples, or products of the other three matrices in the 
plaquettes containing Uj^, and arises from the differentiation of the Wilson action 
S g . Of course, the staples at the border of the lattice contain links denned by 
the boundary conditions (C* or periodic, depending on the direction). The final 
solution for Hj^ is 







Uj^Vj^ - -II- hl , 



(3.22) 



TA 



where the subscript TA indicates the traceless anti-Hermitian part of the matrix: 

(3.23) 



Q TA = i(Q-O t )-itr(Q-gt) . 



3.2. Reducing the number of flavours: the Hybrid algorithm 

Because of Eq. (2.8), the HMC algorithm that we have just described simulates 
eight fermion flavours in the continuum. In order to come back to four flavours, wc 
note that 

Pf(A) =±(dctA) 1 / 2 (3.24) 

(we saw in Eq. (2.35) that det^4 is a positive number). 

One can simulate (det A) 1 / 2 by reverting to an approximate algorithm. A suit- 
able choice is the R algorithm, 5 in which discretization errors in the molecular 
dynamics part are of 0(At 2 ). 

On a lattice closed with C*-BC, the full QCD action in the presence of Nf 
families of degenerate continuum fermions can be written as follows: 

Z = j{VU) [det(AiA) e ] N '' 6 e- s ^ = J (VU) e -s g {u)+^-tr\o^ A) e ^ (3 25) 

where (A^A) e is the restriction of A^A to the even lattice sites. 

Since the implementation of the R algorithm in the present case reduces to 
simple adaptation of a standard technique to the system described by the equation 
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of motion obtained in the previous subsection, we do not elaborate further on this 
point. However, we still have the "sign problem" of Eq. (3.24). The sign of the 
Pfaffian can be included by reweighting the expectation values according to 

_ (0-signPf(A)) + 
{U} - (signPf(A)) + ' (3 - 26) 

where (■■■)+ means that the expectation values have been obtained simulating 
+(det A) 1 / 2 . One then needs to monitorize the sign of the Pfaffian depending on 
the gauge configurations. Some techniques to do that are explained in Ref. 9. 

4. Conclusions 

C* boundary conditions are interesting to study some spontaneous symmetry break- 
ing aspects of QCD. They are relevant when one analyses confinement through 
monopole condensation. We have shown in this work how these boundary condi- 
tions can be implemented to carry out a lattice simulation of full QCD with stag- 
gered fermions. We have proved that the common even-odd trick used to avoid the 
fermion redoubling produced by the introduction of the pseudofermionic field can be 
applied also to this case. However, there is an additional redoubling which forces to 
work with a minimum number of eight flavours with the usual Hybrid Monte Carlo 
algorithm, which we have adapted to this case. An alternative to avoid that is to 
consider a non-exact algorithm, which can be applied to any number of flavours. 

These algorithms have been implemented and are presently running on an APE 
Quadrics machine to explore monopole condensation in full QCD. 
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Appendix A 

We will show that detA is a real number, or, equivalently, that detA* = detA. 
This property is required in order that one can use the matrix A^A to introduce 
the pseudofermionic field. 

Let us consider A in the form 

. ( B —to — D T \ .. 
A= {m + D St J. ( A -!) 

where Dij = Mij — mdij. 

Using the general property of the determinant of a square block matrix reported 
in Eq. (2.30) and exchanging''' the two rows of blocks of A we can rewrite 

= dct(m + D)dct(-m- D T - B(D + m)- 1 B jf ) . (A.2) 
Remembering that = —D, we can write A* in the form 

*={- m % V). 

so that, exchanging the two columns of blocks in A^ and using again Eq. (2.30), we 
can write 

det A T = det „ , nT 

y B —m + D 1 J 

= dct(m- D)dct(-m + D T - B(-D + m)- 1 B*>) . (A.4) 

After extracting a factor (— l)^™ 2 ^ = m 2N from both detA and detA' 1 ' and 
defining a = 1/m, we can write 

detA = m 2N det(l + aD)det(l + aD T +a 2 B(l + aD)- 1 B*) 
det ^ = m 2N det{l-aD)det{l-aD T +a 2 B{l-aD)~ 1 B^) (A.5) 

It clearly appears from Eq. (A.5) that det A* is obtained from det A by changing 
the sign of a. Therefore, in order to show that detA = detA*, it is sufficient to 
show that det A is an even function of a. 

^This corresponds to an even number of row exchanges, since A has dimension 2N and TV, being 
the dimension of the fermion matrix with periodic boundary conditions, is always an even number. 
Therefore det A does not change under this operation. 
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Let us consider dct(l + aD) at first. We can expand the determinant as follows: 
det(l + a£>) = cxp (trln(l + aD)) 

= exp^-trf^t^(aD)^ . (A.6) 

ft is easy to see that the trace of the product of an odd number of D matrices is 
zero. Indeed D only connects nearest neighbour lattice sites, so it is not possible to 
connect a site to itself using the product of an odd number of D matrices. Therefore 
only even powers of a appear in the expansion in Eq. (A.6) and det(l + aD) is an 
even function of a. 

Let us consider now det(l + aD T + a 2 B(l + aD)' 1 B^), which we rewrite as 
det(l + P(a)), where 

P(a) = aD T + a 2 B ( f^(-l) k a k D k J B f . (A.7) 

\fc=0 / 

We notice that the matrix P{a) is expressed as series expansion in a, where the 
coefficient of the fc-th term is a homogeneous polynomial of degree k in the matrices 
B, B^ ,D and D T . Therefore, expanding again the determinant as 

det(l + P(a)) = exp(trln(l + P(a))) 

= exp^-trf;^P(a)^, (A.8) 

we see that det(l + P(a)) can be expanded as a power series in a and that the 
coefficient of the fc-th term is the trace of a homogeneous polynomial of degree k in 
the matrices B,B\D and D T : since these matrices only connect nearest neighbour 
sites, the trace is zero for k odd. Therefore also in this case the determinant is an 
even function of a. 

We conclude that det A, being the product of even functions of a, is also an even 
function of a, and therefore, from Eq. (A. 5), det A = det At. 

Appendix B 

We will give here a proof, alternative to that presented in Ref. 8, of the equality 
of the determinants of the submatrices on even and odd sites of M^M in the case 
with standard boundary conditions. 

In the standard case the fermion matrix M, defined in Eq. (2.3), has the fol- 
lowing form 

M= (|l *m~)' (R1) 
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and it is easily shown that D\ Q = —D oc . Using this last property it follows that 

M.M=(- 2 -^ rf-jLO- <B - 2) 

Using the decomposition of the determinant given in Eq. (2.30), it is easy to show 
that 

det M = det (m 2 - ^D eo D oe ^j . (B.3) 

On the other hand, if we exchange even and odd variables in M before applying 
Eq. (2.30), we obtain 

det M = det (m 2 - ^D oc D co ^j , (B.4) 



and therefore 



det ( m 2 - -D oc D co j = det (m 2 - ^D co D oc ) . ( B.o : 



